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Abstract 

Observable consequences of the hypothesis that the observed universe is a numerical simulation 
performed on a cubic space-time lattice or grid are explored. The simulation scenario is first 
motivated by extrapolating current trends in computational resource requirements for lattice QCD 
into the future. Using the historical development of lattice gauge theory technology as a guide, 
we assume that our universe is an early numerical simulation with unimproved Wilson fermion 
discretization and investigate potentially-observable consequences. Among the observables that are 
considered are the muon g — 2 and the current differences between determinations of a, but the 
most stringent bound on the inverse lattice spacing of the universe, 10 11 GeV, is derived 

from the high-energy cut off of the cosmic ray spectrum. The numerical simulation scenario could 
reveal itself in the distributions of the highest energy cosmic rays exhibiting a degree of rotational 
symmetry breaking that reflects the structure of the underlying lattice. 
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I. INTRODUCTION 



Extrapolations to the distant futurity of trends in the growth of high-performance com- 
puting (HPC) have led philosophers to question — in a logically compelling way — whether 
the universe that we currently inhabit is a numerical simulation performed by our distant 
descendants [lj. With the current developments in HPC and in algorithms it is now pos- 
sible to simulate Quantum Chromodynamics (QCD), the fundamental force in nature that 
gives rise to the strong nuclear force among protons and neutrons, and to nuclei and their 
interactions. These simulations are currently performed in femto-sized universes where the 
space-time continuum is replaced by a lattice, whose spatial and temporal sizes are of the 
order of several femto-meters or fermis (1 fm = 10~ 15 m), and whose lattice spacings (dis- 
cretization or pixelation) are fractions of fermis Q This endeavor, generically referred to as 
lattice gauge theory, or more specifically lattice QCD, is currently leading to new insights 
into the nature of matter [^J Within the next decade, with the anticipated deployment of 
exascale computing resources, it is expected that the nuclear forces will be determined from 
QCD, refining and extending their current determinations from experiment, enabling pre- 
dictions for processes in extreme environments, or of exotic forms of matter, not accessible 
to laboratory experiments. Given the significant resources invested in determining the quan- 
tum fluctuations of the fundamental fields which permeate our universe, and in calculating 
nuclei from first principles (for recent works, see Refs. [JH6]), it stands to reason that future 
simulation efforts will continue to extend to ever-smaller pixelations and ever-larger vol- 
umes of space-time, from the femto-scale to the atomic scale, and ultimately to macroscopic 
scales. If there are sufficient HPC resources available, then future scientists will likely make 
the effort to perform complete simulations of molecules, cells, humans and even beyond. 
Therefore, there is a sense in which lattice QCD may be viewed as the nascent science of 
universe simulation, and, as will be argued in the next paragraph, very basic extrapolation 
of current lattice QCD resource trends into the future suggest that experimental searches 
for evidence that our universe is, in fact, a simulation are both interesting and logical. 

There is an extensive literature which explores various aspects of our universe as a simula- 
tion, from philosophical discussions [lj, to considerations of the limits of computation within 
our own universe [7J, to the inclusion of gravity and the standard model of particle physics 
into a quantum computation [8], and to the notion of our universe as a cellular automa- 
ton [9H12]. There have also been extensive connections made between fundamental aspects of 
computation and physics, for example, the translation of the Church- Turing principle [TBlfTl] 
into the language of physicists by Deutsch |15]. Finally, the observational consequences due 
to limitations in accuracy or flaws in a simulation have been considered |16) . In this work, we 
take a pedestrian approach to the possibility that our universe is a simulation, by assuming 
that a classical computer (i.e. the classical limit of a quantum computer) is used to simulate 
the quantum universe (and its classical limit), as is done today on a very small scale, and ask 
if there are any signatures of this scenario that might be experimentally detectable. Further, 
we do not consider the implications of, and constraints upon, the underlying information, 
and its movement, that are required to perform such extensive simulations. It is the case 
that the method of simulation, the algorithms, and the hardware that are used in future 
simulations are unknown, but it is conceivable that some of the ingredients used in present 

1 Surprisingly, while QCD and the electromagnetic force are currently being calculated on the lattice, the 
difficulties in simulating the weak nuclear force and gravity on a lattice have so far proved insurmountable. 

2 See Refs. [5J |3] for recent reviews of the progress in using lattice gauge theory to calculate the properties 
of matter. 
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FIG. 1. Left panel: linear fit to the logarithm of the CRRs of MILC asqtad and the SPECTRUM 
anisotropic lattice ensemble generations. Right panel: extrapolation of the fit curves into the future, 
as discussed in the text. The blue (red) horizontal line corresponds to lattice sizes of one micron 
(meter), and vertical bands show the corresponding extrapolated years beyond 1999 for the lattice 
generation programs. 



day simulations of quantum fields remain in use, or are used in other universes, and so 
we focus on one aspect only: the possibility that the simulations of the future employ an 
underlying cubic lattice structure. 

In contrast with Moore's law, which is a statement about the exponential growth of 
raw computing power in time, it is interesting to consider the historical growth of mea- 
sures of the computational resource requirements (CRRs) of lattice QCD calculations, and 
extrapolations of this trend to the future. In order to do so, we consider two lattice gener- 
ation programs: the MILC asqtad program jTT] , which over a twelve year span generated 
ensembles of lattice QCD gauge configurations, using the Kogut-Susskind [18J (staggered) 
discretization of the quark fields, with lattice spacings, b, ranging from 0.18 to 0.045 fm, and 
lattice sizes (spatial extents), L, ranging from 2.5 to 5.8 fm, and the on-going anisotropic 
program carried out by the SPECTRUM collaboration j 19 j . using the clover- Wilson [20j |2T] 
discretization of the quark fields, which has generated lattice ensembles at b ~ 0.1 fm, with 
L ranging from 2.4 to 4.9 fm [22] • At fixed quark masses, the CRR of a lattice ensem- 
ble generation (in units of petaFLOP-years) scales roughly as the dimensionless number 
^qcdL 5 /b 6 , where Xqcd = 1 fm is a typical QCD distance scale. In fig. [I] (left panel), the 
CRRs are presented on a logarithmic scale, where year one corresponds to 1999, when MILC 
initiated its asqtad program of 2 + 1 -flavor ensemble generation. The bands are linear fits 
to the data. While the CRR curves in some sense track Moore's law, they are more than 
a statement about increasing FLOPS. Since lattice QCD simulations include the quantum 
fluctuations of the vacuum and the effects of the strong nuclear force, the CRR curve is a 
statement about simulating universes with realistic fundamental forces. The extrapolations 
of the CRR trends into the future are shown in the right panel of fig. [T] The blue (red) 
horizontal line corresponds to a lattice of the size of a micro-meter (meter), a typical length 
scale of a cell (human), and at a lattice spacing of 0.1 fm. There are, of course, many caveats 
to this extrapolation. Foremost among them is the assumption that an effective Moore's 
Law will continue into the future, which requires technological and algorithmic developments 
to continue as they have for the past 40 years. Related to this is the possible existence of 
the technological singularity [231 [21], which could alter the curve in unpredictable ways. 
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And, of course, human extinction would terminate the exponential growth [lj. However, 
barring such discontinuities in the curve, these estimates are likely to be conservative as 
they correspond to full simulations with the fundamental forces of nature. With finite re- 
sources at their disposal, our descendants will likely make use of effective theory methods, 
as is done today, to simulate every-increasing complexity, by, for instance, using meshes 
that adapt to the relevant physical length scales, or by using fluid dynamics to determine 
the behavior of fluids, which are constrained to rigorously reproduce the fundamental laws 
of nature. Nevertheless, one should keep in mind that the CRR curve is based on lattice 
QCD ensemble generation and therefore is indicative of the ability to simulate the quantum 
fluctuations associated with the fundamental forces of nature at a given lattice spacing and 
size. The cost to perform the measurements that would have to be done in the background 
of these fluctuations in order to simulate — for instance — a cell could, in principle, lie on a 
significantly steeper curve. 

We should comment on the simulation scenario in the context of ongoing attempts to 
discover the theoretical structure that underlies the Standard Model of particle physics, and 
the expectation of the unification of the forces of nature at very short distances. There 
has not been much interest in the notion of an underlying lattice structure of space-time 
for several reasons. Primary among them is that in Minkowski space, a non-vanishing spa- 
tial lattice spacing generically breaks space-time symmetries in such a way that there are 
dimension-four Lorentz breaking operators in the Standard Model, requiring a large num- 
ber of fine-tunings to restore Lorentz invariance to experimentally verified levels [2H]- The 
fear is that even though Lorentz violating dimension four operators can be tuned away at 
tree-level, radiative corrections will induce them back at the quantum level as is discussed 
in Refs. [261 [27]. This is not an issue if one assumes the simulation scenario for the same 
reason that it is not an issue when one performs a lattice QCD calculation The under- 
lying space-time symmetries respected by the lattice action will necessarily be preserved at 
the quantum level. In addition, the notion of a simulated universe is sharply at odds with 
the reductionist prejudice in particle physics which suggests the unification of forces with a 
simple and beautiful predictive mathematical description at very short distances. However, 
the discovery of the string landscape [28, 29], and the current inability of string theory to 
provide a useful predictive framework which would post-diet the fundamental parameters of 
the Standard Model, provides the simulators (future string theorists?) with a purpose: to 
systematically explore the landscape of vacua through numerical simulation. If it is indeed 
the case that the fundamental equations of nature allow on the order of 10 500 solutions [30J, 
then perhaps the most profound quest that can be undertaken by a sentient being is the ex- 
ploration of the landscape through universe simulation. In some weak sense, this exploration 
is already underway with current investigations of a class of confining beyond-the-Standard- 
Model (BSM) theories, where there is only minimal experimental guidance at present (for 
one recent example, see Ref. |31j). Finally, one may be tempted to view lattice gauge the- 
ory as a primitive numerical tool, and that the simulator should be expected to have more 
efficient ways of simulating reality. However, one should keep in mind that the only known 
way to define QCD as a consistent quantum field theory is in the context of lattice QCD, 
which suggests a fundamental role for the lattice formulation of gauge theory. 

Physicists, in contrast with philosophers, are interested in determining observable con- 
sequences of the hypothesis that we are a simulation 00 In lattice QCD, space-time is 

3 Current lattice QCD simulations are performed in Euclidean space, where the underlying hyper-cubic 
symmetry protects Lorentz invariance breaking in dimension four operators. However, Hamiltonian lattice 
formulations, which are currently too costly to be practical, are also possible. 

4 There are a number of peculiar observations that .could be attributed to our universe being a simulation, 



replaced by a finite hyper-cubic grid of points over which the fields are defined, and the 
(now) finite-dimensional quantum mechanical path integral is evaluated. The grid breaks 
Lorentz symmetry (and hence rotational symmetry), and its effects have been defined within 
the context of a low-energy effective field theory (EFT), the Symanzik action, when the lat- 
tice spacing is small compared with any physical length scales in the problem [55| U4J [^J 
The lattice action can be modified to systematically improve calculations of observables, by 
adding irrelevant operators with coefficients that can be determined nonperturbatively. For 
instance, the Wilson action can be (9(6)-improved by including the Sheikholeslami-Wohlert 
term [21]. Given this low-energy description, we would like to investigate the hypothesis 
that we are a simulation with the assumption that the development of simulations of the 
universe in some sense parallels the development of lattice QCD calculations. That is, early 
simulations use the computationally "cheapest" discretizations with no improvement. In par- 
ticular, we will assume that the simulation of our universe is done on a hyper-cubic grid H 
and, as a starting point, we will assume that the simulator is using an unimproved Wilson 
action, that produces 0(b) artifacts of the form of the Sheikholeslami-Wohlert operator in 
the low-energy theory 

In section [TTJ, the simple scenario of an unimproved Wilson action is introduced. In 



section III, by looking at the rotationally- invariant dimension- five operator arising from 
this action, the bounds on the lattice spacing are extracted from the current experimental 
determinations, and theoretical calculations, of g — 2 of the electron and muon, and from 
the fine-structure constant, a, determined by the Rydberg constant. Section IV considers 
the simplest effects of Lorentz symmetry breaking operators that first appear at 0(b 2 ), 
and modifications to the energy- momentum relation. Constraints on the energy-momentum 
relation due to cosmic ray events are found to provide the most stringent bound on b. We 
conclude in section IVl 



II. UNIMPROVED WILSON SIMULATION OF THE UNIVERSE 

The simplest gauge invariant action of fermions which does not contain doublers is the 
Wilson action, 



S (w) = b 4 j2 t {w \x) = hi { m+ i) v> WW 



+ y E [ " X ) U ^ ^ x + W ~ + 1) Ul{x - bfi) ^(x -&£)], (1) 

X 

which describes a fermion, ip, of mass m interacting with a gauge field, A^x), through the 
gauge link, 

/ rx+bfi \ 

U^(x) = exp I ig / dzA^(z) , (2) 



but that cannot be tested at present. For instance, it could be that the observed non- vanishing value of 
the cosmological constant is simply a rounding error resulting from the number zero being entered into a 
simulation program with insufficient precision. 
5 Hsu and Zee [32] have suggested that the CMB provides an opportunity for a potential creator /simulator 
of our universe to communicate with the created/simulated without further intervention in the evolution 
of the universe. If, in fact, it is determined that observables in our universe are consistent with those that 
would result from a numerical simulation, then the Hsu-Zee scenario becomes a more likely possibility. 
Further, it would then become interesting to consider the possibility of communicating with the simulator, 
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where jl is a unit vector in the /i-direction, and g is the coupling constant of the theory. 
Expanding the Lagrangian density, £S W \ in the lattice spacing (that is small compared with 
the physical length scales), and performing a field redefinition |3H], it can be shown that the 
Lagrangian density takes the form 

£W = ^pifxp + 77#V + Cp^a^G^ + 0(b 2 ), (3) 

where G^ u = —i[D^ , D u ]/g is the field strength tensor and is the covariant derivative, m 
is a redefined mass which contains 0(b) lattice spacing artifacts (that can be tuned away). 
The coefficient of the Pauli term ijja^G^ip is fixed at tree level, C p = 1 + 0(a), where 
a = g 2 /(4n). It is worth noting that as is usual in lattice QCD calculations, the lattice 
action can be 0{b) improved by adding a term of the form 8C^ W ^ = C sw^ipo pLU G^ 1 ' ip to the 
Lagrangian with C sw = —C p + 0(a). This is the so-called Sheikholeslami-Wohlert term. 
Of course there is no reason to assume that the simulator had to have performed such an 
improvement in simulating the universe. 

III. ROTATION ALLY INVARIANT MODIFICATIONS 

Lorentz symmetry is recovered in lattice calculations as the lattice spacing vanishes when 
compared with the scales of the system. It is useful to consider contributions to observables 
from a non-zero lattice spacing that are Lorentz invariant and consequently rotationally 
invariant, and those that are not. While the former type of modifications could arise from 
many different BSM scenarios, the latter, particularly modifications that exhibit cubic sym- 
metry, would be suggestive of a structure consistent with an underlying discretization of 
space-time. 

1. QED Fine Structure Constant and the Anomalous Magnetic Moment 

For our present purposes, we will assume that Quantum Electrodynamics (QED) is simulated 
with this unimproved action, eq. Q. The 0(b) contribution to the lattice action induces an 
additional contribution to the fermion magnetic moments. Specifically, the Lagrange density 
that describes electromagnetic interactions is given by eq. ([3]), where the interaction with 
an external magnetic field B is described through the covariant derivative = 8^ + ieQA^ 
with e > and the electromagnetic charge operator Q, and where the vector potential 
satisfies V x A = B. The interaction Hamiltonian density in Minkowski-space is given by 

H in t = ^pA^ifS* - id^Qip + ^cr^F^ + C^Ua F^-ty + ... . (4) 
2m 4m 4 

where F^ v = d il A v — d v A ii is the electromagnetic field strength tensor, and the ellipses denote 

terms suppressed by additional powers of b. By performing a non-relativistic reduction, the 

or even more interestingly, manipulating or controlling the simulation itself. 

6 The finite volume of the hyper-cubic grid also breaks Lorentz symmetry. A recent analysis of the CMB 
suggests that universe has a compact topology consistent with two compactified spatial dimensions and 
with a greater than 4tr deviation from three uncompactified spatial dimensions [35J . 

7 The concept of the universe consisting of fields defined on nodes, and interactions propagating along the 
links between the nodes, separated by distances of order the Planck length, has been considered previously, 
e.g. see Ref. (36]. 

8 It has been recently pointed out that the domain-wall formulation of lattice fermions provides a mecha- 
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last two terms in eq. [4] give rise to Hint, mag = —f^ B, where the electron magnetic moment 
H is given by 

H = ^(g + 2mbC p + ...) S = 9(b) ^S, (5) 

where g is the usual fermion g-factor and S is its spin. Note that the lattice spacing 
contribution to the magnetic moment is enhanced relative to the Dirac contribution by one 
power of the particle mass. 

For the electron, the effective g-factor has an expansion at finite lattice spacing of 

+ Ohadrons + 0-u,t + ^weak + f^eb C p + ... , (6) 

where the coefficients Ci, in general, depend upon the ratio of lepton masses. The calculation 
by Schwinger provides the leading coefficient of C2 = \. The experimental value of gi^pt/2 = 
1.001 159 652 180 73(28) gives rise to the best determination of the fine structure constant 
a (at b = 0) [39] . However, when the lattice spacing is non-zero, the extracted value of a 
becomes a function of b, 

a(b) = a(0) - 2nm e b C p + O (a 2 b) , (7) 

where a(0) _1 = 137.035 999 084(51) is determined from the experimental value of electron 
g-factor as quoted above. With one experimental constraint and two parameters to deter- 
mine, a and b, unique values for these quantities cannot be established, and an orthogonal 
constraint is required. One can look at the muon g — 2 which has a similar QED expansion 
to that of the electron, including the contribution from the non-zero lattice spacing, 

+ 4aLns + + «weik + ^ C p + ... . (8) 

Inserting the electron g — 2 (at finite lattice spacing) gives 

^ = ^® + K " m e )b C p + O (a 2 b) . (9) 

Given that the standard model calculation of g^(0) is consistent with the experimental 
value, with a ~ 3.6cr deviation, one can put a limit on b from the difference and uncer- 
tainty in theoretical and experimental values of g£L/2 = 1.001 165 920 89(54)(33) and 

^0^/2 = 1.001 165 918 02(2) (42) (26) [39]. Attributing this difference to a finite lattice 
spacing, these values give rise to 

r 1 = (3.6 ±1.1) x 10 7 GeV , (10) 

which provides an approximate upper bound on the lattice spacing. 



nism by which the number of generations of fundamental particles is tied to the form of the dispersion 
relation |37| . Space-time would then be a topological insulator. 
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2. The Rydberg Constant and a 

Another limit can be placed on the lattice spacing from differences between the value of a 
extracted from the electron g — 2 and from the Rydberg constant, R^. The latter extraction, 
as discussed in Ref . [39J , is rather complicated, with the value of the Roo obtained from a x 2 - 
minimization fit involving the experimentally determined energy-level splittings. However, to 
recover the constraints on the Dirac energy-eigenvalues (which then lead to Roo), theoretical 
corrections must be first removed from the experimental values. To begin with, one can 
obtain an estimate for the limit on b by considering the differences between a's obtained 
from various methods assuming that the only contributions are from QED and the lattice 
spacing. Given that it is the reduced mass (/i ~ m e ) that will compensate the lattice spacing 
in these QED determinations (for an atom at rest in the lattice frame), one can write 

5a = 2iim e b C p , (11) 

where C p is a number 0(1) by naive dimensional analysis, and is a combination of the con- 
tributions from the two independent extractions of a. There is no reason to expect complete 
cancellation between the contributions from two different extractions. In fact, it is straight- 
forward to show that the 0(b) contribution to the value of a determined from the Rydberg 
constant is suppressed by a 4 ml, and therefore the above assumption is robust. In addition to 
the electron g — 2 determination of fine structure constant as quoted above, the next precise 
determination of a comes form the atomic recoil experiments, or 1 = 137.035 999 049(90) q] 
|39j . given an a priori determined value of the Rydberg constant. This gives rise to a differ- 
ence of |5a| = (1.86 ± 5.51) x 10~ 12 between two extractions, which translates into 

b = | (-0.6 ± 1.7) x 10~ 9 | GeV" 1 . (12) 

As this result is consistent with zero, the la values of the lattice spacing give rise to a limit 
of 

r x > 4 x 10 8 GeV , (13) 

which is seen to be an order of magnitude more precise than that arising from the muon 
9-2. 

For more sophisticated simulations in which chiral symmetry is preserved by the lattice 
discretization, the coefficient C p will vanish or will be exponentially small. As a result, 
the bound on the lattice spacing derived from the muon g — 2 and from the differences 
between determinations of a will be significantly weaker. In these analyses, we have worked 
with QED only, and have not included the full electroweak interactions as chiral gauge 
theories have not yet been successfully latticized. Consequently, these constraints are to be 
considered estimates only, and a more complete analysis needs to be performed when chiral 
gauge theories can be simulated. 

IV. ROTATIONAL SYMMETRY BREAKING 

While there are more conventional scenarios for BSM physics that generate deviations in 
g— 2 from the standard model prediction, or differences between independent determinations 



9 Extracted from a 87 Rb recoil experiment |40j . 
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of a, the breaking of rotational symmetry would be a solid indicator of an underlying space- 
time grid, although not the only one. As has been extensively discussed in the literature, 
another scenario that gives rise to rotational invariance violation involves the introduction of 
an external background with a preferred direction. Such a preferred direction can be defined 
via a fixed vector, |41j . The effective low-energy Lagrangian of such a theory contains 
Lorentz covariant higher dimension operators with a coupling to this background vector, and 
breaks both parity and Lorentz invariance [12]. Dimension three, four and five operators, 
however, are shown to be severely constrained by experiment, and such contributions in the 
low-energy action (up to dimension five) have been ruled out [25l W\\ H3J HI] . 

3. Atomic Level Splittings 

At 0(b 2 ) in the lattice spacing expansion of the Wilson action, that is relevant to describing 
low-energy processes, there is a rotational-symmetry breaking operator that is consistent 
with the lattice hyper-cubic symmetry, 

C RV = C™«±^ j^D.D, V , (14) 

where the tree- level value of C RV = 1. In taking matrix elements of this operator in the 
Hydrogen atom, where the binding energy is suppressed by a factor of a compared with 
the typical momentum, the dominant contribution is from the spatial components. As each 
spatial momentum scales as m e a, in the non-relativistic limit, shifts in the energy levels are 
expected to be of order 

5E ~ C RV a A m 3 e b 2 . (15) 

To understand the size of energy splittings, a lattice spacing of b^ 1 = 10 8 GeV gives an 
energy shift of order 5E ~ 10~ 26 eV, including for the splittings between substates in 
given irreducible representations of SO(3) with angular momentum J > 2. This magnitude 
of energy shifts and splittings is presently unobservable. Given present technology, and 
constraints imposed on the lattice spacing by other observables, we conclude that there is 
little chance to see such an effect in the atomic spectrum. 

4- The Energy-Momentum Relation and Cosmic Rays 

Constraints on Lorentz- violating perturbations to the standard model of electroweak interac- 
tions from differences in the maximal attainable velocity (MAV) of particles (e.g. Ref. [25]). 
and on interactions with a non-zero vector field (e.g. Ref. [45]). have been determined pre- 
viously. Assuming that each particle satisfies an energy-momentum relation of the form 
Ef = \pi\ 2 c 2 + m 2 cf (along with the conservation of both energy and momentum in any 
given process), if c 7 exceeds c e ±, the process 7 — > e + e~ becomes possible for photons with 

an energy greater than the critical energy i? C rit. = 2m e c 2 c 7 / \J c 2 ~ c 2 ± , and the observa- 
tion of high energy primary cosmic photons with -E 7 < 20 TeV translates into the constraint 
c 7 — c e ±< 10~ 15 . Ref. [23] presents a series of exceedingly tight constraints on differences be- 
tween the speed of light between different particles, with typical sizes of 5cij< 1CT 21 — 1CT 22 
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FIG. 2. The energy surface of a massless, non-interacting Wilson fermion with r = 1 as a function 
of momentum in the x and y directions, bounded by —tt < bp X)V < ir, for p z = is shown in blue. 
The continuum dispersion relation is shown as the red surface. 



for particles of species i and j. At first glance, these constraints [26] would appear to 
also provide tight constraints on the size of the lattice spacing used in a simulation of the 
universe. However, this is not the case. As the speed of light for each particle in the dis- 
cretized space-time depends on its three-momentum, the constraints obtained by Coleman 
and Glashow |25| do not directly apply to processes occurring in a lattice simulation. 

The dispersion relations satisfied by bosons and Wilson fermions in a lattice simulation 
(in Minkowski space) are 
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respectively, where r is the coefficient of the Wilson term, Eb and Ef are the energy of a 
boson and fermion with momentum k, respectively. The summations are performed over 
the components along the lattice Cartesian axes corresponding to the x,y, and z spatial 
directions. The implications of these dispersion relations for neutrino oscillations along one 
of the lattice axes have been considered in Ref. [46] . Further, they have been considered as 
a possible explanation [47J of the (now retracted) OPERA result suggesting superluminal 
neutrinos |3H]- The violation of Lorentz invariance resulting from these dispersion relations 
is due to the fact that they have only cubic symmetry and not full rotational symmetry, as 
shown in fig. [2j It is in the limit of small momentum, compared to the inverse lattice spacing, 
that the dispersion relations exhibit rotational invariance. While for the fundamental parti- 



cles, the dispersion relations in eq. (16) and eq. (17) are valid, for composite particles, such 



as the proton or pion, the dispersion relations will be dynamically generated. In the present 
analysis we assume that the dispersion relations for all particles take the form of those in 
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eq. (16) and eq. (17). It is also interesting to note that the polarizations of the massless 
vector fields are not exactly perpendicular to their direction of propagation for some direc- 
tions of propagation with respect to the lattice axes, with longitudinal components present 
for non-zero lattice spacings. 

Consider the process p — > p + 7, which is forbidden in the vacuum by energy-momentum 
conservation in special relativity when the speed of light of the proton and photon are equal, 
c p = Cry. Such a process can proceed in- medium when v p > c 7 , corresponding to Cerenkov 
radiation. In the situation where the proton and photon have different MAV's, the absence 
of this process in vacuum requires that \c p — c 7 |< 10~ 23 [25} H9]. In lattice simulations of 
the universe, this process could proceed in the vacuum if there are final state momenta 
which satisfy energy conservation for an initial state proton with energy E,i moving in some 
direction with respect to the underlying cubic lattice. Numerically, we find that there are 
no final states that satisfy this condition, and therefore this process is forbidden for all 



proton momentum _ In contrast, the process 7 — > e + e~ , which provides tight constraints 
on differences between MAV's [25], can proceed for very high energy photons (those with 
energies comparable to the inverse lattice spacing) near the edges of the Brillouin zone. 
Further, very high energy 7r°'s are stable against tt° — > 77, as is the related process 7 — > 7r°7. 

With the dispersion relation of special relativity, the structure of the cosmic ray spec- 
trum is greatly impacted by the inelastic collisions of nucleons with the cosmic microwave 
background (CMB) [5TH 1ST] . Processes such as 7cmb + N — > A give rise to the predicted 
GKZ-cut off scale [5"U|I5"T] of ~ 6x 10 20 eV in the spectrum of high energy cosmic rays. Recent 
experimental observations show a decline in the fluxes starting around this value [52, 53J, 
indicating that the GKZ-cut off (or some other cut off mechanism) is present in the cosmic 
ray flux. For lattice spacings corresponding to an energy scale comparable to the GKZ cut 
off, the cosmic ray spectrum will exhibit significant deviations from isotropy, revealing the 
cubic structure of the lattice. However, for lattice spacings much smaller than the GKZ 
cut off scale, the GKZ mechanism cuts off the spectrum, effectively hiding the underlying 
lattice structure. When the lattice rest frame coincides with the CMB rest frame, head-on 
interactions between a high energy proton with momentum |p| and a photon of (very- low) 
energy oj can proceed through the A resonance when 

^ 2 |P| 2 (^, n ^ 



id 



4|p| 



1 + 
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- m \-< br + ... , (18) 
4|p| v ' 

for |p| <C 1/b, where 9 and are the polar and azimuthal angles of the particle momenta in 
the rest frame of the lattice, respectively. This represents a lower bound for the energy of 
photons participating in such a process with arbitrary collision angles. 

The lattice spacing itself introduces a cut off to the cosmic ray spectrum. For both the 
fermions and the bosons, the cut off from the dispersion relation is E max ~ 1/b. Equating 
this to the GKZ cut off corresponds to a lattice spacing of b ~ 1CT 12 fm, or a mass scale of 
b^ 1 ~ 10 11 GeV. Therefore, the lattice spacing used in the lattice simulation of the universe 
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A more complete treatment of this process involves using the parton distributions of the proton to relate 

its energy to its momentum [26J. For the composite proton, the p — > p + 7 process becomes kinematically 

allowed, but with a rate that is suppressed by C(Aq CD 6 7 ) due to the momentum transfer involved, 

effectively preventing the process from occuring. With momentum transfers of the scale ~ 1/b, the final 

states that would be preferred in inclusive decays, p — > + are kinematically forbidden, with invariant 

masses of ~ 1/b. More refined explorations of this and other processes are required. 
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must be 10~ 12 fm in order for the GZK cut off to be present or for the lattice spacing 
itself to provide the cut off in the cosmic ray spectrum. The most striking feature of the 
scenario in which the lattice provides the cut off to the cosmic ray spectrum is that the 
angular distribution of the highest energy components would exhibit cubic symmetry in the 
rest frame of the lattice, deviating significantly from isotropy. For smaller lattice spacings, 
the cubic distribution would be less significant, and the GKZ mechanism would increasingly 
dominate the high energy structure. It may be the case that more advanced simulations will 
be performed with non-cubic lattices. The results obtained for cubic lattices indicate that 
the symmetries of the non-cubic lattices should be imprinted, at some level, on the high 
energy cosmic ray spectrum. 

V. CONCLUSIONS 

In this work, we have taken seriously the possibility that our universe is a numerical simula- 
tion. In particular, we have explored a number of observables that may reveal the underlying 
structure of a simulation performed with a rigid hyper-cubic space-time grid. This is mo- 
tivated by the progress in performing lattice QCD calculations involving the fundamental 
fields and interactions of nature in femto-sized volumes of space-time, and by the simulation 
hypothesis of Bostrom [T]. A number of elements required for a simulation of our universe 
directly from the fundamental laws of physics have not yet been established, and we have 
assumed that they will, in fact, be developed at some point in the future; two important 
elements being an algorithm for simulating chiral gauge theories, and quantum gravity. It 
is interesting to note that in the simulation scenario, the fundamental energy scale defined 
by the lattice spacing can be orders of magnitude smaller than the Planck scale, in which 
case the conflict between quantum mechanics and gravity should be absent. 

The spectrum of the highest energy cosmic rays provides the most stringent constraint 
that we have found on the lattice spacing of a universe simulation, but precision measure- 
ments, particularly the muon g — 2, are within a few orders of magnitude of being sensitive 
to the chiral symmetry breaking aspects of a simulation employing the unimproved Wilson 
lattice action. Given the ease with which current lattice QCD simulations incorporate im- 
provement or employ discretizations that preserve chiral symmetry, it seems unlikely that 
any but the very earliest universe simulations would be unimproved with respect to the 
lattice spacing. Of course, improvement in this context masks much of our ability to probe 
the possibility that our universe is a simulation, and we have seen that, with the excep- 
tion of the modifications to the dispersion relation and the associated maximum values of 
energy and momentum, even 0(b 2 ) operators in the Symanzik action easily avoid obvious 
experimental probes. Nevertheless, assuming that the universe is finite and therefore the 
resources of potential simulators are finite, then a volume containing a simulation will be 
finite and a lattice spacing must be non-zero, and therefore in principle there always remains 
the possibility for the simulated to discover the simulators. 
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